function [dQ, n, q] = calpathVar(alphar,alphap,rad,V)
global R0 Vs g0 S m
alt = (rad-1)*R0;
[~, sonic, ~, rho] = atmoscoesa(alt);
Ma = V*Vs/sonic;
[CL, ~] = getCLCD(alphar, Ma);
[~, CD] = getCLCD(alphap, Ma);
q = 0.5*rho*(V*Vs)^2;
L = q*S*CL/m;
D = q*S*CD/m;
n = sqrt(L^2 + D^2)/g0;
dQ = 6E-8*sqrt(rho)*(V*Vs)^3.15;
end